Four-body continuum-discretized coupled-channels calculations using a transformed 

harmonic oscillator basis 
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The scattering of a weakly bound three-body system by a target is discussed. A transformed 
harmonic oscillator basis is used to provide an appropriate discrete and finite basis for treating 
the continuum part of the spectrum of the projectile. The continuum-discretized coupled-channels 
framework is used for the scattering calculations. The formalism is applied to different reactions, 
6 He+ 12 C at 229.8 MeV, 6 He+ 64 Zn at 10 and 13.6 MeV, and 6 He+ 208 Pb at 22 MeV, induced by the 
Borromean nucleus 6 He. Both the Coulomb and nuclear interactions with a target are taken into 
account. 
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I. INTRODUCTION 

The development of radioactive nuclear beam facili- 
ties has allowed the study of nuclei far from the line 
of stability, bringing to the fore new nuclear structure 
problems. A significant topic in recent years has been 
the study of halo nuclei [J, H, 0|. These are weakly 
bound, spatially extended systems, typically comprising 
a core and one or two valence nucleons. Particularly in- 
teresting examples of such systems are Borromean nu- 
clei, i.e., three-body composite systems with no binary 
bound states. These nuclei have attracted special atten- 
tion because their loosely bound nature reflects a delicate 
interplay between two- and three-body forces, constitut- 
ing a challenge to existing theories, and a motivation for 
the development of new ones. The detailed structure of 
the continuum spectrum of these systems is still not fully 
understood, partially due to the ambiguities associated 
with the underlying forces between the constituents. Due 
to their low binding energy, halo nuclei are easily broken 
up in the nuclear and Coulomb field of the target nu- 
cleus. Therefore few-body reaction theories, developed 
to extract reliable information from experimental data 
of reactions involving loosely bound systems, have to in- 
clude, as an essential ingredient, a realistic description of 
coupling to the continuum part of the spectrum. 

From the theoretical point of view, the treatment of re- 
actions involving loosely bound systems must deal with 
the complication that these continuum breakup states 
are not square-normalizable. A convenient method to 
circumvent this problem is to replace the states in the 
continuum by a finite set of normalized states, thus pro- 
viding a discrete basis that, hopefully, can be truncated 
to a small number of states and yet provide a reliable 
description of the continuum. Several prescriptions to 
construct such a discrete basis have been proposed. For 



two-body composite systems, where the continuum states 
are easily calculated, one can use a discretization proce- 
dure in which the continuum spectrum is truncated at a 
maximum excitation energy and divided into energy in- 
tervals. For each interval, or bin, a normalizable state is 
constructed by superposition of scattering states within 
that bin interval. The method, normally used in the 
continuum-discretized coupled-channels (CDCC) frame- 
work has been very useful in the description of elas- 
tic and breakup observables in reactions involving weakly 
bound two-body projectiles. 

An alternative to the binning procedure is to represent 
the continuum spectrum by the eigenstates of the internal 
Hamiltonian in a basis of square integrable (or L 2 ) func- 
tions, such as Laguerre 0,0, @], Gaussian or Stur- 
mian (ill. fl^. fl3l. [l4| functions. In practice, the diagonal- 
ization is performed in a finite (truncated) set of states 
and the resulting eigenstates, also known as pseudo-states 
(PS) , are regarded as a finite and discrete representation 
of the spectrum of the system. The pseudo-states are 
then used within a coupled-channels calculation in the 
same way as the continuum bins. 

The PS method has the appealing feature of being 
readily applicable also to describe the spectrum of three- 
body systems, in which case the Hamiltonian is diagonal- 
ized in a complete set of square-integrable functions for 
the three-body Hilbert space. Several applications of this 
method can be found in the literature, for both structure 
[TBI ] and reaction problems [10]. In the latter case, the 
method is an extension of the CDCC formalism to re- 
actions with three-body projectiles, using a pseudo-state 
model for the continuum. 

One such PS method proposed recently is the Trans- 
formed Harmonic Oscillator (THO) method [U H]. 
Given the ground-state wave function of the system, the 
THO method performs a Local Scale Transformation 
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(LST) [19] that converts the bound ground-state wave 
function of the system into the ground-state wave func- 
tion of a Harmonic Oscillator (HO). Once the LST is 
obtained, the HO basis can be transformed, by the in- 
verse LST, to a discrete basis in the physical space. The 
THO basis functions are not eigenfunctions of the Hamil- 
tonian (except for the ground state) but the Hamiltonian 
can be diagonalized in an appropriate truncated basis 
to produce approximate eigenvalues and eigenfunctions. 
This method has been shown to be useful for describing 
the two-body continuum in both structure [l?], [13, Hi 
and scattering [2ll H3, HUl problems. In a recent work 
[HI the THO method was generalized to describe contin- 
uum states of three-body systems, based on expansion 
in Hyperspherical Harmonics (HH) jlH]. In particular 
the method was applied to the Borromean nucleus 6 He, 
for which several strength functions, including the dipole 
and quadrupole Coulomb transition strengths, were cal- 
culated. These observables are found to converge quickly 
with respect to the number of THO basis states included. 
Furthermore, the calculated strength distributions are in 
very good agreement with those obtained using three- 
body scattering wave functions [2^ |. 

Most of our knowledge of 6 He comes from the analy- 
sis of reactions where secondary beams collide with sta- 
ble nuclei. These experiments have been performed with 
both light [13, Hi and heavy targets (H, M ED, E3, El, 
and at low and high energies, providing a body of data 
which can be used to benchmark reaction and structure 
models. The theoretical understanding of reactions in- 
volving a three-body projectile, such as 6 He, is a compli- 
cated task because it requires the solution of a four-body 
scattering problem. At high energies, a variety of approx- 
imations have been used such as semiclassical approxima- 
tions [34l. l35l. l36| . frozen halo or adiabati c ap proximations 
[37l.l38l|. Multiple Scattering expansions [33, H, 52] , four- 
body DWBA [42l |43|. among others. However, at ener- 
gies of a few MeV per nucleon, some of these approxima- 
tions are not justified. Then the use of the CDCC method 
is an alternative to solve these problems. For a four-body 
problem (three-body projectile) this method has already 
been applied using a PS basis based on Gaussian func- 
tions. The scattering of 6 He by 12 C [lj| and 209 Bi Q 
have been studied. In both cases a good agreement was 
obtained with the experimental data of Refs. [H,|4^| and 
[3lj . respectively. 

In this work, we study the scattering of a three-body 
projectile by a target using the CDCC formalism. The 
novel feature of the present approach is the use of the 
THO PS basis to represent the states of the projectile. 
These states are then used to generate the projectile- 
target coupling potentials that enter the system of cou- 
pled equations. Furthermore, we have developed a new 
procedure to calculate these coupling potentials making 
use of an expansion of the wave functions of the projectile 
internal states in a HH basis. 

This paper is structured as follows. In Section II the 
three-body discretization method is presented. In Sec- 



tion III the multipole expansion of the interaction poten- 
tial between the projectile and the target is addressed. 
In Section IV we describe the three-body model for the 
Borromean nucleus 6 He. In Section V we apply the for- 
malism to the reactions 6 He+ 12 C at £i ab =229.8 MeV, 
6 He+ 64 Zn at £j ab =13.6 and 10 MeV, and 6 He+ 208 Pb at 
Slab =22 MeV. Finally, Section VI summarizes and draws 
conclusions. 



II. THREE-BODY CONTINUUM 
DISCRETIZATION METHOD 

The THO discretization method applied to a three- 
body system is described in detail in Ref. [2^ . For com- 
pleteness, in this Section we outline the main features of 
the formalism. In the three-body case, it is convenient 
to work with the hyperspherical coordinates {p,a,x,y}. 
They are obtained from the Jacobi coordinates {x, y} 
that are illustrated in Fig. [TJ The variable x is propor- 
tional to the relative coordinate between two of the parti- 
cles, with a scaling factor depending on their masses [IS] 
and y is proportional to the coordinate from the center 
of mass of these two particles to the third particle, again 
with a scaling factor depending on their masses. From 
these coordinates, the hyperradius (p) and the hyperan- 
gle (a) are defined as p — y 1 x 1 + y 2 and tana = x/y. 
Obviously there are three different Jacobi sets but p is 
the same for all of them. 

For a three-body system the discretization method has 
two parts. First, the wave functions of the system are 
expanded in Hyperspherical Harmonics (HH) [25| . We 
define states of good total angular momentum as 

X ^ lm M^ V ) T KL l Ws x , (1) 

where T^-^(O) are the hyperspherical harmonics that 
depend on the angular variables £1 = {a,x, y}, Xs^ ls 
the spin wave function of the two particles related by 
the coordinate x, and x L i ls the spin function of the 
third particle. Each component of the wavefunction 
(or channel) is defined by the set of quantum numbers 
(3 = {K,l x ,l v ,l, S x , jab}- Here, K is the hypermomen- 
tum, l x and l y are the orbital angular momenta associ- 
ated with the Jacobi coordinates x and y, I = l x + l y is 
the total orbital angular momentum, S x is the spin of the 
particles related by the coordinate x, and j ab = I + S x . 
Finally, j = j ab + I is the total angular momentum, with 
I the spin of the third particle, which we assume fixed. 
The physical states of the system can now be expressed 
as a linear combination of the states given by Eq. {1} as 

il> ili (j>,n) = ^Rfi j (p)yf>j ll (n), (2) 

where {Rpj} are the hyperradial wave functions. 



Secondly, the THO method is used to obtain the func- 
tions Rpj(p). Writing the ground-state wave function in 
the form of Eq. the equation that defines the LST 
for each channel (3 is 
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where R B p{p) is the bound ground-state hyperradial 
wave function for the channel /3, with N B p the normal- 
ization factor, and Rqx ( s ) is the ground-state hyperra- 
dial wave function of the HO for the hypermomentum 
K, that is already normalized. Finally, the THO basis is 
constructed for each channel by applying the LST, sp(p), 
to the HO basis 
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where the L^(t) are generalized Laguerre polynomials 
and M%K is the normalization constant of the HO basis. 
Here the index i denotes the number of hyperradial exci- 
tations. Note that as i increases, the functions Rfp (p) 
become more oscillatory and explore larger distances. 

For channels with quantum numbers that do not con- 
tribute to the ground-state wave function, the (ground 
state) channel with the closest quantum labels to the 
channel of interest is used to construct the LST. One im- 
portant point concerns the label K which governs the p K 
behavior of the hyperradial wave function close to the 
origin. To guarantee the correct behavior of the wave- 
function, we select a channel from the ground-state wave 
function with the same K . If this is not possible, a chan- 
nel with K — 1 is used and the corresponding hyperradial 
wave function is then multiplied by p. 

The required discrete eigenstates are now calculated by 
diagonalizing the three-body Hamiltonian of the projec- 
tile in a finite THO basis up to rib hyperradial excitations 
in each channel, 



b THO 
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where n labels the eigenstates for a given angular momen- 
tum j and e n j will be the associated energy. Replacing in 
this expression the functions i)™? {p, 0) by their explicit 
expansion in terms of the HH, Eq. |(5]), and performing 
the sum in the index i for i = 0, . . . , rib, we can express 
the PS basis states as 



(7) 







Note that the choice of the HO parameter has no in- 
fluence in the calculation of the LST since changes to 
this parameter are equivalent to making a linear trans- 
formation in the oscillator variable s. This gives the same 
result for the right part of Eq. <j3j) . 
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Figure 1: (Color online) Relevant coordinates for the scatter- 
ing of a three-body projectile by a structureless target. 



III. MULTIPOLE EXPANSION OF THE 
PROJECTILE- TARGET POTENTIAL 

The eigenstates given in Eq. |(7|) are a discrete repre- 
sentation of the states of the three-body projectile. From 
them, the four-body wavefunction of the projectile-target 
system, schematically depicted in Fig. [U is formed as 
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where R is the coordinate from the target to the center of 
mass of the projectile, L is the orbital angular momentum 
of the projectile-target relative motion and J is the total 
angular momentum, J = L + j. The radial functions 
finj (R) satisfy the system of coupled equations 
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fLj( R ) 
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(9) 



where m r is the reduced mass of the projectile-target 
system. The coupling potentials V^ n j L , n ,j, (R) are then 

VL 3 ,L'n>A R ) = {LnjJM\V pt {Ti,r 2 ,r3)\L'n'j'JM), 

(10) 

where the ket \LnjJM) denotes the function 

®Lni( R > x >y) s iven b y 



$ J L %(R,x,y) 



J2 4>n" (^ y)(LM L] p\JM)Y LML (R). 

(11) 

To calculate these coupling potentials, a multipole ex- 
pansion of the projectile-target interaction is developed. 
The procedure is analogous to that for a three-body 
problem reported in Ref. In that work the tradi- 

tional method of bin averaging was used as discretization 
method instead of the THO method. We assume that the 
projectile-target interaction is the sum of the interactions 
of each particle of the projectile with the target, VktiXk) 
with k = 1,2,3. For each pair potential, an appropriate 
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Jacobi set is chosen so that the corresponding coordinate 
r k depends only on the vectors R and y k . Assuming that 
the potentials are central, the coefficients of the multipole 
expansion are generated as 



K(R,y k ) 



1 



+1 



V k (r k )P Q (z k )dz k , 



(12) 



where Pg(z k ) is a Legendre polynomial, Q is the mul- 
tipole order and z k = y k R is the cosine of the angle 
between y k and R. So, the coupling potential can be 
expressed as 



V£ njfL , n , f (R)=Y,(- 1 ) J ~ iiL '( L o)w(LLijf,QJ)F^ nlj ,(R) 

r 



Q 



(13) 



where the radial form factor F9 , -,(R) is 
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(sin a k ) 2 (cos a k ) 2 da k p 5 dp Rl™ (p)<p l ^ k (a k )V h Q {R, y k )^f k (a k )R^P f (p), (14) 
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with /3fc being the set of quantum numbers in the fc'th 
Jacobi system where the potential depends on x k , and 
/? being the set in the Jacobi system in which the states 
of the projectile are calculated. The matrix elements 
Npf3 k transform the hyperangular, angular and spin part 
of the wave functions from one Jacobi set to another. 
Their explicit expression as a function of the Raynal- 
Revai coefficients is developed in Ref. Note that 

Eqs. lfl~3|) and (fl"4"l) are completely general, and do not 
depend on the nature of the basis. 



IV. STRUCTURE MODEL FOR HE 

The 6 He nucleus is treated here as a three-body sys- 
tem, comprising an inert a core and two valence neu- 
trons. The ground state has total angular momentum 
jir _ q+ w ^h experimental binding energy of 0.973 MeV. 
The ground state wave function was obtained by solv- 
ing the Schrodinger equation in hyperspherical coordi- 
nates, following the procedure described in [2 a . I48II. and 
making use of the codes FACE [H]+ STURMXX |49l|. In 
these calculations, the n- 4 He potential was taken from 
Ref. It consists of an energy independent Woods- 

Saxon potential, supplemented by a spin-orbit term with 
a Woods-Saxon derivative radial shape. This potential 



reproduces the low-energy s- and p-phase shifts up to 10 
MeV. For the NN interaction we used the potential pro- 
posed by Gogny, Pires and Tourreil (GPT) [5l|, which 
contains central, spin-orbit and tensor components. This 
interaction was developed to give simultaneously an ac- 
ceptable fit to two nucleon scattering data up to 300 MeV 
and to describe reasonably the properties for finite nu- 
clei, particularly the radii, within the Hartree-Fock ap- 
proximation. Besides the two-body (n — n and n — a) 
potentials, the model Hamiltonian also includes a simple 
phenomenological three-body force, depending only on 
the hyperradius, according to the following power form 



V 3 b(p) = - 



i + (p/by 



(15) 



where a, 6, and c are adjustable parameters. This po- 
tential is introduced to correct the under-binding caused 
by our neglect of other configurations, such as the t+t 
channel. 

We have performed different calculations that truncate 
the maximum hypermomentum at K max — 2,4,6,8, and 
10, respectively. For each value of K max , the three-body 
potential has been adjusted to give the same binding en- 
ergy and mean square radius (for j* = + states) and the 
same position for the 2 + resonance (for = 2+ states). 
The latter value was also used for the j* = 1~ states. 
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Figure 2: Energy spectrum of the states with j n = + , = 
1 , and 7 = up to 30 MeV excitation, obtained for a 
THO basis with A" m ax = 8 and rib = 4. 



The parameter a varies with i^ max and j, being of the 
order of 4 MeV for f = 0+ and 3 MeV for f = 1~ , 2+, 
The parameter b varies with K max , within the range 4-6 
fm. The parameter c was fixed to 3 in all cases. 

The number of channels (3 for each calculation in- 
creases drastically with K max , making the calculations 
much more demanding computationally. In the follow- 
ing, unless stated otherwise, the calculations presented 
use the basis with K max = 8. As we will show below, 
this basis provides converged results with respect to the 
hypermomentum for all the reactions considered in this 
work. For this case, the number of channels j3 is 15 for 
f = 0+, 26 for f = \~ and 46 for f = 2+. The cal- 
culated three-body wave function has a binding energy 
of 0.95372 MeV and a rms point nucleon matter radius 
of 2.46 fm when assuming an alpha-particle rms matter 
radius of 1.47 fm. 

The Jacobi set in which the two neutrons are related 
by the coordinate x is chosen to generate the THO ba- 
sis. Applying the THO method and diagonalizing the 
Hamiltonian in a finite THO basis, a set of eigenstates 
is obtained. For j = + , the diagonalization produces 
a state with negative energy, that corresponds to the 
ground state of the system. The remaining eigenvalues 
appear at positive energy, and are then associated with 
a discrete representation of the continuum spectrum. As 
an example, in Fig.[2]we present the distribution of eigen- 
values obtained for a basis with rib = 4, for the states 
j = 0+, j = 1~, and j = 2+ up to 30 MeV. 



V. APPLICATION TO REACTIONS 

In this Section we apply the formalism developed in 
Sec. IIII I to different reactions induced by 6 He, taking 
the eigenstates obtained with the THO basis to repre- 
sent the projectile states. We note that, even for a small 
value of rib, the THO method will produce eigenstates 



at very high excitation energies. States above a certain 
excitation energy will not be relevant for the description 
of the collision process, since they will be very weakly 
coupled. For this reason, in these calculations the basis 
is truncated at a maximum excitation energy, and only 
those eigenstates below this value were included in the 
coupled-channels calculation. The maximum energy is 
chosen independently for each reaction and each rib, in 
order to achieve convergence of the results with respect 
to this energy. 

In the present calculations, only the 6 He states with 
j = + ,l~,2 + are considered. Previous CDCC calcula- 
tions for the reactions studied in this work 
have shown that using these partial waves is sufficient 
to obtain converged results and to describe satisfactorily 
the existing data for such reactions. We did not attempt 
to include higher partial waves, since this would make 
the calculations very demanding computationally . 

For these included f states the coupling potentials 
given by Eq. lfl3|) are calculated for multipolarities Q = 
0, 1,2. Both Coulomb and nuclear interactions were in- 
cluded. We emphasize that continuum-continuum cou- 
plings were also included. The diagonal as well as non- 
diagonal coupling potentials were generated by folding 
the neutron-target and cv-target interactions according 
to Eq. (fT0|) . These interactions are represented by phe- 
nomenological optical potentials at the relevant projec- 
tile incident energy per nucleon [53]. Then, the coupled 
equations ([9]) are solved with the code FRESCO (54T |. that 
reads the coupling potentials from external files. In most 
cases, the Numerov method was used to solve the coupled 
equations. However, in some cases, particularly when 
excitation energies close to the total kinetic energy are 
involved, this method was found to be numerically un- 
stable, and the R-matrix method [H3| was used instead. 
This method is more time consuming but has the advan- 
tage of being numerically more stable. In the following, 
we present the results for different reactions for which 
experimental data exist. 

6 He+ 12 C. We study this reaction at 229.8 MeV, for 
comparison with the experimental data of Lapoux et 
al. j45|. The n+ 12 C potential was taken from the global 
parametrization of Watson et al. [56]. The a+ 12 C poten- 
tial was represented in terms of a standard Woods-Saxon 
shape with the parameters adjusted in order to reproduce 
the elastic data for this system at 34.75 MeV per nucleon 
The parameters for these potentials are listed in 
Table H 

The coupled equations were solved up to J = 70 and 
the solutions were matched to their asymptotic form at 
the radius R m = 200 fm. In Fig. [3] we present the an- 
gular distribution of the elastic differential cross section 
relative to Rutherford. The thick solid line is the full 
CDCC result for a basis with nb=4. This calculation re- 
produces the data fairly well (open circles) up to 10°, 
but it clearly underestimates the data points at larger 
angles. Interestingly, this effect was also found in the 
phenomenological analysis of Lapoux et al. [45j], as well 
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Table I: Optical model parameters used in this work. All po- 
tentials are parametrized using the usual Woods-Saxon form, 
with a real volume part and volume (W v ) and surface (Wd) 
imaginary part. Reduced radii are related to physical radii 

by R = r A 1 T / ' i . 



System 


V 


ro 


Q() 


Wv 


w d 


n 


a; 




(MeV) 


(fm) 


(fm) 


(MeV) 


(MeV) 


(fm) 


(fm) 


n+ 12 C 


49.46 


1.115 


0.57 


3.05 


7.48 


1.15 


0.5 


a+ 12 C 


100. 


1.289 


0.71 


19.98 




1.738 


0.495 


n+ 64 Zn 


51.82 


1.203 


0.668 


0.29 




1.203 


0.668 












5.97 


1.279 


0.534 


a+ 64 Zn 


123 


1.676 


0.43 


20.40 




1.467 


0.43 


n+ 208 Pb 


47.37 


1.222 


0.726 




6.24 


1.302 


0.351 


a+ 208 Pb 


96.44 


1.376 


0.625 




32. 


1.216 


0.42 



as in the four-body CDCC calculation of Matsumoto et 
al. [lB| for the same reaction. We also show the analo- 
gous calculation when omitting all the couplings to the 
continuum (one channel calculation) with a dashed line. 
For the reaction at 229.8 MeV we conclude that the effect 
of coupling to the continuum is a reduction of the cross 
section for angles beyond 5°. This effect has also been 
observed in the scattering of 11 Be+ 12 C at E ~ 49 MeV 
per nucleon [32} , and is probably present in other reac- 
tions induced by weakly bound projectiles at energies of 
a few tens of MeV per nucleon. That the no-continuum 
coupling calculation reproduces the data reasonably well 
at the larger angles is probably fortuitous, and cannot be 
attributed to the adequacy of this approximation. As we 
have shown, continuum couplings are very important in 
this reaction. 

We also show in Fig. [3] the full CDCC calculation for 
«b=2 (dotted line). This calculation is practically undis- 
tinguishable from the calculation with Ub=4, indicating 
that it is not necessary in this case to have a very pre- 
cise discretization of the continuum in terms of excitation 
energy. We found that a maximum excitation energy of 
£max =30 MeV provided good convergence for all the val- 
ues of rib presented. 

6 He+ 64 Zn. We have studied this reaction at two dif- 
ferent energies, namely, 13.6 MeV and 10 MeV, for which 
experimental data exist [58]. The n+ 64 Zn potential was 
taken from the global parametrization of Koning and De- 
laroche [H§|. For the a+ 64 Zn system, we took the opti- 
cal potential derived in Ref. [58j]. The parameters are 
listed in Table [1 The coupled equations were solved up 
to ,/ = 60 and 40, respectively, and for projectile-target 
separations up to R m = 100 fm. 

In Figs. [4] and [5] we present the experimental and calcu- 
lated angular distributions of the elastic cross section for 
these two reactions. The dashed lines correspond to the 
one channel calculations (i.e., omitting the continuum) 
and the thick solid lines are the full four-body CDCC 
calculations for a basis with rib =4. 

At £i ab =13.6 MeV (Fig. 0}, the one-channel calcula- 
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Figure 3: (Color online) Elastic differential cross section rela- 
tive to Rutherford as a function of the scattering angle in the 
projectile-target center of mass for the reaction 6 He+ 12 C at 
iJi a b=229.8 MeV. Experimental data are from Ref. [451 ] . 



tion exhibits a pronounced rainbow peak at around 30°, 
which is much smaller in the data. Also, this calculation 
gives a too small cross section at large angles. Inclu- 
sion of couplings to the continuum suppresses this rain- 
bow, and enhances the backward angles cross section, 
improving the agreement with the data in the whole an- 
gular range. In the same figure, we also show the full 
CDCC calculation for a basis with rib=2 (dotted) and 
6 (dot-dashed). These two calculations are very close 
to rifc=4 showing a very good convergence with respect 
rib. The maximum excitation energy required for con- 
vergence depended somewhat on the value of nb, ranging 
from e max =7 MeV (for nb=2) to e max =6 MeV (for nb=6) 

At £iab=10 MeV (Fig.©, the full CDCC calculation 
also improves the agreement with the data at backward 
angles, although some underestimation remains. Inter- 
estingly, the data suggests the presence of a rainbow at 
around 50°, which is not present in our calculation. It 
should be noted that the experimental error bars are large 
at this energy, so more accurate measurements would be 
needed to make more definite conclusions about this ap- 
parent discrepancy. Again, in the same figure, we show 
the full CDCC calculation for a basis with U(,=2 (dotted) 
and 6 (dot-dashed). In this case, we find that the conver- 
gence with respect rib is slower. However the calculations 
with rib=4 and 6 are quite close and give a reasonable 
convergence. As in the previous case, the maximum ex- 
citation energy required for convergence depended some- 
what on the value of nb, ranging from £ max =9 MeV (for 
nb=2) to e max =5 MeV (for nb=6). 

6 He+ 20S Pb. We have performed calculations for this 
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Figure 4: (Color online) Elastic differential cross section rel- 
ative to Rutherford as a function of the cm. scattering angle 
for the reaction 6 He+ 64 Zn at £i a b=13.6 MeV. Experimental 
data are from Ref. [5Sl |. 
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Figure 5: (Color online) Elastic differential cross section rel- 
ative to Rutherford as a function of the cm. scattering angle 
for the reaction 6 He+ 64 Zn at -Ei a b=10 MeV. Experimental 
data are from Ref. [58|]. 



reaction at 22 MeV, in order to compare with the recent 
data of Sanchez-Bemtez et al. We took the n+ 208 Pb 
potential from Ref. ffnj and the a+ 208 Pb potential from 
Ref. [62|. The parameters for these potentials are also 
listed in Table HI The coupled equations were solved up 
to J = 150 and matched to their asymptotic solution at 




30 60 90 120 150 180 



Figure 6: Convergence of the differential elastic cross section 
with respect to .Kmax, for the reaction 6 He+ 208 Pb at £i a b=22 
MeV. All the calculations use rib = 4 for the number of hy- 
perradial excitations and the maximum excitation energy was 
set to 8 MeV. 



Rm = 200 fm. 

First, we discuss the convergence of the calculation 
with respect the hypermomentum (i^ max ) and the hy- 
perradial excitation (rift). In Fig. [6l we show the calcu- 
lations with different values of -Kmax — 2,4,6,8 and for 
the same value of n&=4. For a meaningful comparison, 
in all these cases the three-body potential was adjusted 
in order to give the same binding energy and rms radius, 
for j = + and the same position for the resonance, for 
j = 2 + . We found a relatively fast convergence with re- 
spect to this parameter. In particular, the calculations 
with .Kmax = 6,8 and 10 are very similar (for clarity, 
the latter has been omitted from the figure). For rest 
of reactions the results are quite similar, achieving the 
convergence for -K max =6 or 8. 

The convergence with respect to rib for this reaction 
is illustrated in Fig. [71 For clarity, we show only the re- 
sults for even values of n&. Unlike the previous cases, 
the convergence rate found in this case was rather slow. 
Although the differences in the calculated cross sections 
are less than 5%, the oscillatory pattern at the rainbow 
region changes from one value of rib to another. A pos- 
sible explanation for this slow convergence rate is given 
below. 

In Fig. [8] we compare the experimental and calculated 
angular distributions of the elastic cross section. The 
dashed line is the one channel calculation and the thick 
solid line the full CDCC calculation including the contin- 
uum. The latter uses K max = 8, e max = 8 MeV, n b = 4. 
The one channel calculation shows a rainbow that disap- 
pears in the full calculation, in agreement with the data. 
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Figure 7: Convergence of the differential elastic cross section 
with respect to m, for the reaction 6 He+ 208 Pb at -Ei a b=22 
MeV. 

At backward angles, the agreement with the data is im- 
proved when we include the coupling to the continuum. 
In order to show the contribution of the couplings to each 
j, we also include in this figure the calculation including 
only j = + states (dotted line) and the calculation with 
j = + ,l~ states (thin solid line). From these calcula- 
tions we can conclude that dipole couplings are the main 
responsible for the characteristic reduction of the cross 
section at the angles around the rainbow. The strong 
influence of dipole couplings might explain the slow con- 
vergence with respect to the parameter rib found for this 
reaction. These couplings are very sensitive to the exci- 
tation energy of dipole states, which appear at different 
positions in our discrete representation of the 6 He contin- 
uum, as we vary the number of hyperradial excitations, 
rib. By contrast, in the 6 He+ 12 C case, dipole excitations 
are very small, and this might explain the fast conver- 
gence with respect to rib in that case. 

Moreover, we find that the range of the form factors 
[Eq. (jT4j)] changes significantly for the different pseudo- 
states as rib is changed. This could also contribute to 
the slow convergence at scattering energies close to the 
Coulomb barrier. 



VI. SUMMARY AND CONCLUSIONS 

The collision of a loosely bound three-body projectile 
with a target nucleus has been studied in the framework 
of the continuum-discretized coupled-channels (CDCC) 
method. A set of normalizable states, also known as 
pseudo-states, is used to represent the three-body contin- 
uum of the projectile. In particular we took the Trans- 




Figure 8: (Color online) Elastic differential cross section rela- 
tive to Rutherford as a function of the scattering angle in the 
projectile-target center of mass for the reaction 6 He+ 208 Pb 
at £iab=22 MeV. The full CDCC calculations uses rib =4 and 
£max = 8 MeV. Experimental data are from Ref. [601 ]. 



formed Harmonic Oscillator (THO) basis, which is con- 
structed from the ground state of the system. Within 
the spirit of the CDCC approach, a multipole expansion 
of the coupling potentials has been developed for a four- 
body system (three-body projectile plus a target). 

The formalism has been applied to the reactions 
6 He+ 12 C at 229.8 MeV, 6 He+ 64 Zn at 13.6 and 10 MeV, 
and 6 He+ 208 Pb at 22 MeV, taking into account both the 
Coulomb and nuclear interactions. 

Overall, we find good agreement between the calcu- 
lated and the experimental elastic scattering angular dis- 
tributions. However, for the 6 He+ 12 C reaction at 229.8 
MeV the calculations underestimate the experimental 
data for cm. scattering angles beyond 10°. The fact 
that this effect was also found in previous analyses of 
this reaction fl6l . |45| suggests that the discrepancy is not 
related to the particular features of our approach. 

For the reaction 6 He+ 64 Zn at 13.6 and 10 MeV the 
calculations are in fair agreement with the data, the re- 
production being better in the higher energy case. At 10 
MeV our calculations do not predict a rainbow at around 
50°, a hint of which is seen in the data, but is broadly 
consistent with the data within the stated experimental 
errors. 

In actual coupled-channels calculations, the discrete 
basis has to be truncated in the excitation energy (s max ), 
the maximum hypermomentum (-Kmax), and the maxi- 
mum number of hyperradial excitations (rib). In all the 
cases under study, we have found a good convergence of 
the calculated observables with respect to the parame- 
ters £ max and -Kn iax . However, the rate of convergence 
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with respect to rib was found to depend very much on 
the specific reaction. For the reaction 6 He+ 12 C at 229.8 
MeV the convergence was found to be very fast, with 
rib = 2 providing fully converged results. For 6 He+ Zn 
at near-barrier energies, we required rib ~ 4 for an ac- 
ceptable convergence. Finally, for 6 He+ 208 Pb at 22 MeV, 
the convergence was found to be slow and oscillatory. In 
fact, our biggest calculation, corresponding to rib = 6, 
is still not fully converged. Because of computational 
limitations we have not explored this question further, 
as required to study the convergence of the calculations 
with respect to the basis size. 

This work shows that the use of the transformed har- 
monic oscillator basis, developed in previous works, com- 
bined with the standard CDCC method, provides a re- 
liable procedure for the treatment of the scattering of a 
loosely bound three-body projectile by a target. It will 
be interesting to compare this method with other repre- 
sentations of the continuum, including the standard dis- 
cretization procedure in terms of continuum bins which, 
in the case of three-body projectiles, requires the calcu- 



lation of the three-body scattering states. This work is 
underway and the results will be published elsewhere. 
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